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ABSTRACT 



A method for analytically determining the stress 
distribution in a moment loaded circular semimonoeoque 
cylinder with cutout is presented. The technique employ- 
ed in the analysis consists of combining the stringers 
and skin into a smaller number of equivalent stringers 
and shear membranes and, after determining the bending "i 

strain energy of the rings \inder various shear patterns, 
deriving the energies of the corresponding idealized 
structures. Analytical results for the necessary patterns 
of three idealizations along with the method of application 
to specific problems are presented in a form suitable 
for machine computation using Castigllano *s Theorem, 

The comparison of the results of application to an exper- 
imental structure with experimental values is also set 
forth. 
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CHAPTER I 



INTRODUCTION, STATEMEOT? OP PROBLEI4, AND ORGANIZATION 

I. INTRODUCTION 

Because of the increasing demand for high performance 
rocket and missile systems, refined structural design and 
stress prediction methods have become very important# 

Most such techniques are, however, complicated and require 
the services of high grade engineering personnel for their 
use# On the other hand, the construction of experimental 
models for tost is equally expensive and, if the results 
lead to redesign, such models are difficult to restructixre# 
Of particular interest are methods which lead to 
reasonable stress prediction for structures with cutouts 
such as are required for instrioment and fuel system access 
ports# Since overall weight is a vital concern, necessary 
reinforcement of these areas must be efficiently done# If 
the stress predictions can be made with sufficient accuracy, 
consistency, and speed, the development time may be sig- 
nificantly reduced# 



II. THE PROBLEM 



It was the purpose of this study (1) to devise a 
standard method of analysis for the preliminary predi<^ion 
of stresses in cylindrical semimonocoque shells with cutout 
when loaded in pure bending; (2) to present the results 
in such a manner that they may be used without excessive 
re-analysis to predict stresses in experimental structures; 
and (3) to compare the resvilts of the analysis with a set 
of experimental results, 

III. ORGANIZATION 



Chapter II will contain a discussion of the previous 
related studies and the differences between them and this 
study* Also included will be a definition of the various 
special terms used. 

Chapter III is a discussion of the structural ideali- 
zation* the general method of analysis * and the steps of 
the structure analysis. 

Chapter IV presents a svunmary of the results of the 
analysis and the author *s conclusions. 

Appendix I illustrates the stepwise analysis of a 
representative structure. 

Appendix II presents the results of a comparison 
between analytical and experimental results. 
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CHAPTER II 



RELATED STUDIES AND DEFINITIONS OF TERMS 
I. RELATED STUDIES 



The method of analysis used In this thesis was sug- 
gested by a method for handling redundant structures 

Introduced by Professor W. R, Jensen in his Advanced 

■* 

Structtares course at Webb Institute of Naval Archltect\are« 
Reference (1) outlines its application to a test cylinder* 
Professor Jensen made his personal notes covering addition- 
al work in the area fvilly available to this author. 

A closely related method of analysis is presented in 
reference (3). Basically it is a pertvu?bation load tech- 
nique in which individual cancellation loads are placed 
at the boundaries of the desired cutout in an analytical 
structure whose external load system corresponds to the 
physical structure. The stress distribution of the struc- 
ture without cutout is presumed to be known and reference 
(4) supplies tables of coefficients to calculate the 
effect of each type of cancellation load throughout the 
structure (i.e. the resulting stress distribution in the 
unloaded structure in the vicinity of the applied unit 
perturbation loads), in application, then, perturbation 
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loads in magnitude and orientation ne.oessary to cancel 
existing stresses in the loaded whole atructvire are ap- 
plied to the botindaries of the desired cutout and their 
associated stress distributions are superimposed on the 
existing pattern. The resultant stresses represent the 
distribution in the loaded cylinder with cutout. Refer- 
ence (5) shows an excellent agreement with the data of 
reference (2). 

In the author’s opinion the two foregoing studies 
are the only ones sufficiently close to the simplified 
analysis being sought to be pertinent. 

II. DEFINITIONS OF SPECIAL TERMS USED 

Semimonocoque . As used herein, this term denotes a 
structure of uniform thickness with regularly spaced 
stiffeners in the form of ring frames and stringers. 

Orthogonal . As used herein, this tex*m designates 
a property of loads or load systems. When more than one 
load is applied to a structure simultaneously the effects 
of each may be considered Independently if they are or- 
thogonal. The energy expression for the loaded structure 
will contain no cross-product terms between the several 
loads or load systems. 
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CHAPTER III 



IDEALIZATION AND METHOD OF ANALYSIS 

I. IDEALIZATION 



The types of structxires with which thl^ thesis is 
concerned (i.e. cylindrical semimonocoque) suggest, by 
their physical appearance, the kind of idealized structure 
which will be used. It will consist of rings, panels, 
and stringers which have the following characteristics: 

1. The rings are loaded by shear flows only which 
act to produce ring bending. 

2. The panels ^e loaded by shear flows only. 

Their thickness is equal to the thickness of 
panels in the actual structure so that they 
will have the same effectiveness for carrying 
the shears as do panels in the actual structure* 

3. The stringers take longitudinal loads and 
longitudinal shear flows and have effective 
cross sectional areas which include the longi- 
tudinal load carrying capacity of the skin of 
the actual structure. 

Three idealizations, differing only in the n\jmber';jOf 
stringers involved, were chosen for analysis. In general# 
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the number of stringers determines the complexity of the 
analysis. For this thesis, idealizations with six, eight, 
and twelve stringers were ehosen. 

II. METHOD OF ANALYSIS 

Having chosen an idealization it was desired to set 
up the minimim number of independent shear flow and force 
patterns which, by their subsequent combinations, would 
produce any possible loading: pattern in the Idealized 
structure . It is to be noted that, once a ring shear 
flow pattern is fixed, the pattern of panel shear flows 
and stringer forces is also fixed. Figure III - 1 illus- 
trates this principle. The problem, then, reduces to 
that of determining the nmber of independent shear flow 
patterns which may be Imposed upon one ring of the ideal- 
ization. Since shears must remain constant on the ring 
segments between stringers, must produce no net vertieal 
or horizontal forces, and must produce no net torque about 
an axis thru the center and perpendlcxilar to the plane of 
the ring, the number of patterns is limited to the total 
number of shear segments minus the number of static equili- 
brium conditions (i.e., in the six stringer ring, 6-5:3 
independent patterns). 

The applied load, pure bending, requires that the 
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Figure III - 1 




Aflaxmed ring shear flow pattern 




Resulting bay shear flow and stringer fore# pattern 



Note that the top and bottom stringers have twice 
the load beeause shear flows act on both sides of 
these stringers* 
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line of intersection of its action plane with the ring 
plane be an axis of symmetry. If this line passes thru 
segments of the ring rather than thru stringer points the 
segments may have no shear. In the six stringer case, 
then, only one pattern of shear flow would be available, 
figure III-2, For this reason, as well as the fact that 
some v/ork (Ref, (1) ) had already been done with symmetry 
thru stringers, it was decided to place the plane of action 
thru stringers. It should be noted that this makes the 
validity of the comparison with the selected experimental 
data (Ref, (2) ) somewhat doubtful. This point will be 
covered in Chapter IV and Appendix II, 

For convenience it was decided to designate individual 
ring shear flow patterns by the type of symmetry possessed 
about the vertical and horizontal axes. With S standing 
for Symmetry and A standing for Antisymmetry, the first 
letter of the hyphenated pair in a parenthesis designates 
the type of symmetry about the vertical axis v;hile the 
second letter designates the type of symmetry about the 
horizontal axis. 

Before choosing the ring shear patterns, one more 
consideration was taken into account - orthogonality of 
load systems. As defined in Chapter II, it is a property 
of orthogonal load systems that the strain energy of a 
structure loaded by several such systems is equal to the 
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Figure III - 2 
SjjmetTj Axis 




This flgxire demonstrates that. If a symmetry axis 
Is put thru the panels, the six stringer Idealization 
is limited to one shear pattern, (S-S), No other 
(S-S) nor (S-A) pattern Is possible which satisfies 
static equllibrlxaa* 
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sum of the strain energies of the structure loaded by 
each separately. The coefficients of terms containing 
load cross products are equal to zero. Patterns with 
different symmetry are orthogonal, but, if a choice existed 
among patterns with the same symmetry, the ones chosen 
were those which were orthogonal. 

In summary, then, ring shear patterns must conform 
to the following conditions: 

1, Static equillhrlura - necessary condition, 

2, Vertical axis of symmetry - necessary condition, 

« 

3, Orthogonality - desireable condition. 

The actual ring shear patterns chosen were arrived 
at by trial and error. An interesting by-product of this 
process was the data for table III-l which gives the 
distribution of ring shear patterns by symmetry type for 
each of the three Idealizations. The proof of this table 
for the six stringer Idealization is given in Appendix I, 
but it was not proven for the eight and twelve stringer 
cases. Table III-l is presented primarily as a point of 
interest. 

Once the ring patterns were chosen, corresponding 
half structures v/ere set up ( the structure is symmetrical 
about its centerline plane ). For a given pattern the 
ring bending strain energy was computed and expressed as 
one-half the square of the shear flow multiplied by a 
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Table III - 1 
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A 


Ring Shear Patterns 


Ideal .Str* 


(S-S) 


(S-A) 


(A-A) 


(A^S) 


6 Stringer 


1 


1 


1 


0 


8 Stringer 


2 


1 


1 


1 


12 Stringer 


3 


2 


2 


2 



l^latrlbutlon of Ring Shear Pattema by Symetry 
Type for Idealized Struettirea* 
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"ITexlbllity factor," This factor is the same for any 
ring loaded with the given pattern. The half- structure 
was considered on a bay by bay basis as demonstrated in 
Appendix I, and the loadings on the successive bays 
starting at the half-bay (the center bay of the overall 
structure) were determined. The strain energy from stringers, 
panels, and rings was determined and the sum gave the 
strain energy of the half-structure. 

If loads are orthogonal, the sum of the strain energies 
for each pattern is equal to the strain energy of the 
idealized structxire loaded by the combined patterns. Or- 
thogonality was checked for patterns of like symmetry. 

III. SUMMARY 



The steps of the analysis are: 

1, Choose an idealization, 

2, Select ring shear patterns. 

3, For each pattern set up the half- structure, 

4, For each pattern calculate the ring "Flexibility 
Factor." 

5, For each pattern, on a bay by bay basis, calcu- 
late the energy matrix factors. 

Note that the combination of the patterns is not done in 
the analysis. This step is taken when the results are 
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applied to a given physical structure since only then 
will the boxmdary conditions be known (i.e. the required 
relationships between the pattern shears and forces and 
the forces applied to the structure). The complete analy- 
sis of a particular idealization consists of the matrices 
which correspond to the ring shear patterns for the ideal- 
ization, These are found in Chapter IV for each of the 
three idealizations. 
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CHAPTER IV 



SUMMARY OP RESULTS AND CONCLUSIONS 
I. SUMI4ARY OP RESULTS 

The nature of the results indicates that the best 
method of presentation is in the form of figures and 
energy matrices* They ape so presented on pages 17 to 
32 . It is intended by this form of presentation to make 
the results available in the most compact and useable 
form possible. The figures which accompany the results 
are intended to assist the reader in visualizing the 
structure to which the matrix applies. The six stringer 
idealization is presented in its entirety for both sym- 

I 

metry patterns while, for the eight and twelve stringer 
structures, only a sample bay (the half-bay, bay mmber 
(1 ) ) is illustrated for each pattern. Note that stringei^ 
force values are diminished as they propogate thru the 
structure and terminate in a force designated by a capital 
letter. The end of the structure is considered as built 
into a rigid wall so that the moment of inertia of the 
end "ring" is considered to approach infinity (i.e. its 
flexibility factor is zero). The meaning of the physical 
constants is summarized in table IV - 1. 



The idealizations of this thesis are all for a nine 
bay structure because the most readily available test 
data (Ref. (2) ) were taken from such a structure. If 
application is desired to a structvire having a different 
number of bays, the matrices must be altered. The step 
by step analysis of one pattern of the six stringer idea- 
lization is presented in Appendix I so that a model is 
available for anyone wishing to make such alterations. 
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TABLE IV - 1 



TABLE OP MATRIX SYMBOL DEFINITIONS 
q - Shear flow. 

R - Radius to the neutral axis of the skin. 

( PR - Radius to the neutral axis of the skin-ring com- 

bination. 

- Length of one bay., 

^ - Thickness of the skin. 

A - Cross-sectional area of the idealized stringer. 

E - Elastic Modulus of the structure material. 

G - Shear Modulus of the structure material. 

^ - Moment of Inertia of a cross-section of the skin- 

ring combination about its neutral axis . 

Underlined terras are taken from the physical struc 

ture. A is computed as follows: 

A=N CAst )+l 2 trRt, 
n n 

where 

N = Number of physical stringers, 
n = Number of idealized stringers, 

Ast = Cross-sectional area of the physical stringer 
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In •quatlon form, th« (S-S) pattern strain enargy Is 
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Each of tha terms In tha above energy expression 
contains a eoafflalent times a second degree load product. 
Such an arrangement suggests tha use of the less "cluttered" 
matrix form and, since modern digital computers are often 
adjusted to solve matrix Inputs, this form was adopted. 
Conversion is accomplished by laying out load rows and 
columns as on page (19C) and placing eaeh coefficient in 
the block corresponding to the intersection of the column 
of one load with the row of the other. Since a matrix is 
symmetric about its major diagonal (i.e. load squared 
terras), non-squared load product positions occur in pairs 
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and one-half of the aoefficient is plaeod In each. For 
the same reason, only one half of the matrix need bo 
filled in since corresponding blocks in the other half 
are Identical. Additional information on general matrix 
\ise may be found in reference (6) while its specif it appll 
cation to structxoral analysis is well Illustrated in refer 
enee (7). 
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II. CONCLUSIONS. 



It was originally hoped that this thesis would pro- 
vide the basis for a standardized method of preliminary 
stress analysis for the type of struoture with which it 
is concerned. Disregarding, for the moment, the compari- 
son plots of Appendix II, the application and advantages 
of the method will first be discussed. 

Even though the simplicity of the Idealization pre- 
cludes great accuracy, this analysis can be used in pre- 
liminary design to provide a relatively quick chock on 
the adequacy of the structure, V/hen the design becomes 
finalized, a detailed chock, using one of the more accur- 
ate methods of stress analysis (e,g, that presented in 
reference (3)) or using an experimental model, may be made. 

In the exact form contained herein, the analysis 
should be applied only to a nine bay structure with load 
axis thru stringers and cutout in the center or center 
plus symmetrically placed adjacent bays, A structural 
analyst would be required to set up the boundary equations 
defined by the physical structure and to derive the elimi- 
nation matrix, but the calculations of matrix factors for 
computer use and transcription of results can be accom- 
plished by less highly trained personnel. The physical 
parameters (l.e, those items underlined on table IV - 1, 
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page (16) may be varied without rederlvlng the elimina- 
tion matrix so that a check on redesign may be rapidly 
made « 

The method of analysis may also be used to obtain 
results corresponding to structures which vary in number 
of bays from that assumed in the thesis. Furthermore, 
by ’’tagging” parameters from a given bay, an analysis 
may be mad# of a structure in which physical parameters 
vary from bay to bay, such as changes in stringer cross 
sections. Finally non-cutout cases may be handled by 
combining the idealizations properly. 

The major objection to the above which may be cited 
is the evidence of poor agreement with experimental results 
seen in Appendix II. In spite of these somewhat disap- 
pointing results, it is felt that the obvious advantages 
warrant a more thorough check than time permitted. It 
would be especially Interesting to have experimental stress 
data in the 180 degree area where the twelve stringer 
idealization predicts a minimum of compression (the rever- 
sal to slight tension actually computed for one of the 
bays is considered highly unlikely). 

In summary, it is felt that, even with its relatively 
poor agreement with the only experimental data available, 
the analysis presented herein could be profitably used in 
preliminary design to quickly obtain approximate structural 
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seantlings. Further, the possible advantages seem to 
the author to warrant an experimental model set up speci- 
fically to verify the analytical data and an Investigation 
of the effect of minor changes in analysis. 
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APPENDIX I 



\ 



DETAn>ED ANALYSIS OF SIX STRINGPJ{ (S-S) PATTERN 



I. AVAILABLE RING PATTERNS 



Plgur# AI - 1 illustratoa the only three possible 
independent ring shear patterns for this idealisation* 
The total number of such patterns is fixed by the three 
static stability conditions at six minus throe equals 
three. An (A~S) pattern is not possible because a hori- 
zontal symmetry axis requires zero shear on the side 
panels while a vertical antisymmetry axis requires a 
continuity of shear across the upper and lower pairs of 
panels. With these conditions it is not possible to 
satisfy both 

ZP„ • 0 

and 

ZM’O. 

Since the load pattern requires a vertical symmetry 
axis, only the (S-S) and (S-A) shear patterns are avail- 
able. 
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II. (S-S) ^STRUCTURE ANALYSIS 



Flgur® AI - 2 and its accompanying calculations ill- 
ustrate the procedure for finding the bay energy as the 
sum of the energies of the elements (i.e. ring, panels, 
and stringers ) . 

Figures AI - 3 thru AI -7 v/lth tables AI - 1 thru 
AI - 5 Illustrate the procedure for finding the success- 
ive bay energies. Table AI - 6 represents the total 
structure energy and is the sum in each matrix block, of 

t 

the corresponding blocks in the matrices of each bay. 
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Figure AI - 3 



BRY © 



F « 







F= ^ 

F = /s “ fscD-^ 

SIX STRINGER (S-S) PATTERN 



7 



( 42 ) 



Table AI -1 
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Figure AI - 5 
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APPENDIX II 



COMPARISOfI WITH }:JCPERIMENTAL RESULTS 



I. COMPARISON STRUCTURE AND PROCEDURE 



Referonoe (3) presents the stress data obtained by 
loading, in pur© bending, an experimental cylindrical 
semlmonocoque shell in which a centrally located cutout 
of varying circumferential length had been made. These 
data wore chosen as the most readily available check on 
the validity of the practical application of this analy- 
sis, Accordingly, the physical constants and combination 
equations necessary to make each analytical idealization 
conform to the experimental structure were derived. The 
combined matrices for each idealization were solved for 
redundante on an IBM 650 computor. Prom the redundant 
solutions, the stringer stress values at several stations 
along the cylinder were computed. For each such station 
a plot was made of analytical and experimental stress 
values versus angular position. The zero angle position 
was vertically downward, figure All - i. 

Because ring and skin shear values in the analytical 
structure are constant between stringei’ positions and 
because the angular extent of such values is so nrach 
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larger than in the experimental structure, it was felt 
that plots of these quantities wovild not be particiolarly 
Informative. An exception to this is the eight stringer 
case where, because the idealization provided a cutout 
size identical to that of an experimental structure, 
panel shear comparisons v/ere made for the bay containing 
the cutout. 

Before the comparison results are presented, three 
qualifying points should be noted: 

1. The axis load symmetry in the experimental 
structure is placed thru the center of panels. 

All analytical results were obtained by placing 
this axis thru stringer points - points of 
natural shear discontinuity. 

2 . In both the six and twelve stringer idealizations, 
the cutout size used was 120 degrees. The 
closest experimental data are for a 130 degree 
cutout . 

3. The moment is so applied to the structure that 
the cutout side is in tension. The same nume- 
rical value of moment is considered to be appli- 
ed to the idealized structure as was applied 

to the physical structure. 
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Figure All - 1 
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II. SIX STRINGER COMPARISON 



The application of a particiilar idealization to a 
moment loaded physical struotvire is accomplished by com- 
bining the various symmetry patterns with the moment 
induced stringer loads so that the loads acting on cut- 
away structure are iliminated. In the six stringer case, 
with a 120 degree cutout, the superposition yields the 
following equations: 

1 , (2 - ;2 - (x-fs - ^ ^S© '^) - F’m - O 

~ =0 

FVn is the moment induced tension load in the bottom string- 
er and equation 1 . is derived from the requirement that 
this stringer have a zero load at each end of the middle 
bay. Equation 2. imposes a zero shear load on the panels 
on either side of the bottom stringer in bay 1. 

The combination equations are used to eliminate 
redundants from the matrix which results from the super- 
position of the pattern matrices. These equations also 
introduce the applied moment in the form of Pm. For the 
six stringer case the final combination equations are: 

1 , s - nq 

2 . = 

Figures All - 2, 3, and 4 give the comparison plots 
with reference data for the 130 degree cutout. 
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III. EIGHT STRINGER COMPARISON 



The three patterns must be superimposed with the 
applied moment to produce zero load in the bottom stringer 
and zero shear in the two adjacent panels. The following 
two elimination equations result, * 

1 , "f S \ ~ ~ ^ 

Figures All - 5, 6, 7, and 8 give the comparison 

■* 

plots with reference data for the 90 degree cutout. 






( 60 ) 




s 



>1 










\ 




( 62 ) 



% 





( 63 ) 




C64) 



IV. TWELVE STRINGER COMPARISON 



In order to achieve a 120 degree cutout, the three 
patterns must be superimposed with the applied moment to 
produce zero load in the bottom stringer and in the two 
lowest side stringers. In addition, the fovir lowest 
shear panels must have zero shear. The following four 
elimination equations result. 

1 * I qS. - 

2* ^5 3<D-C • 

3. -f, 3 = f%i. + fa.^ + 

4* f(xi ~~^si “ ■^S'i " •^o-i 

Figures All - 9, 10, and 11 give the comparison 
plots with reference data for the 130 degree cutout. 
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